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FIRST  LINE  OF  TEXT 


THE  DETERMINATION  OF  THE  ORBIT  OF  AN  ARTIFICIAL 
EARTH  SATELLITE  FROM  THREE  OBSERVATIONS 

0,  M.  Bazhenov 

This  paper  explains  a  method  for  determining  the 
Initial  orbit  of  an  artificial  earth  satellite  from 
three  observati  ons  separated  by  sufficiently  long 
Intervals  of  time. 

To  determine  the  orbit  of  an  AES  (artificial  earth  satellite) 
we  can  use  Gauss »s  well-known  method  which  is  used  to  determine  the 
orbits  of  small  planets  and  comets.  The  successful  application  of 
this  method  is  possible  only  under  conditions  in  which  the  three 
observations  of  a  heavenly  body  are  close  with  respect  to  time,  or 
more  accurately,  under  conditions  in  which  the  change  in  the  mean 
anomaly  in  the  period  of  time  between  observations  will  be  of  the  order 
of  5  to  15°.  In  this  case,  however,  the  determination  of  the  orbit 
is  unreliable  due  to  the  lack  of  precision  in  the  initial  data  obtained 
from  the  observations.  For  those  cases  in  which  the  change  in  the 
mean  anomaly  between  observations  is  great,  the  Gauss  method  becomes 
unwieldy  and  inconvenient . 

The  author  of  the  present  article  published  a  work  in  19J50  [1] 
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(Bazhenov,  1930)  in  which  he  described  a  rather  simple  method  for 

K.  XT 

determining  an  elliptical  orbit  from  three  arbitrary  observations 
which  may  be  separated  by  any  time  intervals  desired,  as  long  as 
they  are  not  too  small. 

This  method  is  satisfactory  in  determining  the  undisturbed 
Keplerian  orbit  of  an  AES  in  cases  where  the  change  in  the  mean 
anomaly  will  be  anywhere  from  5  to  1000°  for  the  period  of  time  be¬ 
tween  observations.  This  method  is  also  used  in  cases  where  the 
variations  in  the  mean  anomaly  are  higher  than  1000°,  however,  in 
these  cases  we  cannot  disregard  the  disturbances  and  assume  that  the 
orbit  is  Keplerian. 

We  present  below  the  above-indicated  method  for  determining  the 
orbit  of  an  AES.  The  observations  of  the  AES  give  the  following 
values : 

t„  «„  8(  (/=  1,  2,  3), 

l.e.  the  moments  of  the  observations,  the  right  ascensions  and  de¬ 
clinations  of  the  AES  corresponding  to  these  moments  and  referring 
to  a  certain  period  (e.g.  1950. 0) .  In  addition,  the  values  9,  X  and 
h  should  be  kncwn,  l.e.  the  geographical  latitude,  longitude  and 
altitude  (in  meters)  above  the  sea  level  of  the  observation  points. 

As  the  unit  of  measurement  for  distances  in  this  work,  we  take 
the  equatorial  radius  of  the  earth  according  to  Krasovskiy  ( aQ  ® 
6,578,245  m)  . 

The  equatorial  rectangular  geocentric  coordinates  of  the  observa¬ 
tion  point  are  determined  from  the  formulas: 


/,  h 

[1  *+*  — 
\ 

|p' cos/  cos  S, 

(»-y 

/  cos  <?'  sin  S, 

(1-3 

p'sin/. 
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p«  and.  <p  are  the  geocentric  radius-vector  and  latitude  of  the  obser- 
vation  point,  while  £3  is  the  local  sidereal  time  (more  accurately, 
the  hour  angle  of  the  point  of  the  vernal  equinox  of  the  period  to 
which  the  coordinate  system  refers) .  The  following  relations  are 
valid: 

p'cos?'  =  1.000339  cos?  —  0.000840  cos  3? 0.000001  cosSo,  \  ,  , 

p' sin?' =  0.995811  sin?  —  0.000336  sin 3? -4- 0.000001  sin5?.  )  '  2' 


The  sidereal  time  may  be  determined  with  a  sufficient  degree 
of  accuracy  from  the  formula 

S=S0  -+-0798565  N-*-  360:98565  n  ^ 

where  SD  is  the  sidereal  time  ( in  degrees)  at  midnight  GMT  for  the 
zero  number  of  the  month  of  the  observation  taken  from  the  table,  N 
is  the  datum  of  the  moment  of  observation  and  n  is  that  fraction  of 
the  twenty-four  hour  periods  from  Greenwich  midnight  to  the  moment 
of  observation. 

Sidereal  time  at  average  midnight  GMT  for  the  zero 
number  of  each  month  In  195?,  1959  and  19b0 


1957 

1958 

1959 

1960 

J  tujuary  .  . .  . . 

99^3940 

9981542 

9S?9142 

9886740 

February  . . 

129.9492 

129.7094 

129.4694 

129.2292 

Karch  . . . . . 

. . 

157.5472 

157  3074 

157.0673 

1578128 

188.1019 

187  8621 

187.6220 

188.2673 

217.6710 

217.4312 

217.1911 

217.9364 

J IUIO  «••••••••••,  a  .  •  .  •  . 

248.2269 

247.9862 

247.7461 

248.4915 

277.7957 

277.5558 

277.3156 

278.0611 

AUgUSt  . 

308.3509 

308.1110 

307 .8708 

338.6162 

C>ep  tac.ber . 

338.9016 

338.6658 

338.4^57 

339.1710 

Oc  tobor  •  •  •  • 

8.2248 

8.2348 

7.9946 

8.7401 

November . . . 

37.0295 

38.7895 

38.5493 

39.2948 

Deceooer . . 

68.5990 

68.3589 

68.1187 

68.6641 

If  we  take  into  consideration  the  shift  in  the  point  of  the  ver¬ 
nal  equinox  along  the  equator  in  the  time  t  =  t  -  1950  ( in  tropical 
years)  which  extends  from  the  moment  of  1950.0  to  the  moment  of 
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observation  t,  then  from  S,  obtained  from  formula  (5),  we  must  subtract 
the  value 

<{•  =  0.0128053  t -4- 0?000000038833  **  (  4) 

and  use  in  the  formulas  ( 1)  the  difference  obtained  instead  of  13. 

A  more  precise  calculation  of  the  precession  can  be  made  in  the 
following  way.  Form  formulas  (1)  we  determine  the  coordinates  X, 

Y  and  Z  in  the  coordinate  system  corresponding  to  the  moment  of 
observation,  and  then  from  the  formula 

(■"^1 650.0 
^1050.0 
650.0/ 

we  obtain  the  coordinates  of  the  observation  point  in  the  system 
corresponding  to  moment  1950. 0.,  The  values  XY,  Y_. . . .  can  be  obtained 
from  Table  5  in  the  work  of  Zagrebin  and  Shumikhina  [5]  (1954)  . 

The  idea  for  the  method  of  calculation  of  sidereal  time  from 
formula  (5)  was  borrowed  from  the  work  of  Zhongolovich,  Amelin  and 
Sabanina  [2]  (1959)* 

The  geocentric  rectangular  equatorial  coordinates  of  the  AES 
are  determined  from  the  formulas: 

x=\?-*-X,  y  —  p-p-4-Y,  *  =  >p -+- Z,  (6) 

In  which  p  Is  the  distance  fran  the  observation  point  to  AES  (the 
unknown  value)  ,  and 

*.  =  cos&cosa,  u.  =  cososin  a,  *  =  iin  S. 

Evidently  x‘  h- p.* -+->*=:  1 . 

In  the  future  we  will  use  of  these  relations! 

4  » 

=  /?*  =  X*  -4-  Y1  -f-  Z*. 

C  =  \X-*-nY-i->Z,  S'  =  R'  —  C', 

p  = _ c  -+-  \Jr 1 _ y*, 


(T) 


(8) 


x ,  — yx  -z; 
-x,  Y,  z, 
-X.  Y.  Z. 


(5) 
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r  and  R  are  the  geocentric  distances  of  the  AES  and  of  the  observation 

r:R..i  LiNt  Of  itJT 

point. 

In  the  method  proposed  In  this  paper  for  determining  the  AES 
orbit,  we  adopt  the  following  as  the  main  unknowns:  a  -  the  major 
simlaxis  of  the  orbit  and  the  two  values: 

g—tanis\ n2,  h  =  —(a*  i  cos  2 ,  (9) 

where  JL  Is  the  angle  of  inclination  of  the  plane  of  the  orbit  of 
the  AES  to  the  plane  of  the  earth* s  equator,  A  is  the  right  ascension 
of  the  ascending  node  of  the  AES  at  the  equator. 

Since  the  undisturbed  AES  orbit  is  planar,  the  following  relation 
is  valid 

xg-+-yh-*~  z  —  0.  (10) 

The  following  equality  follows  from  expression  ( 10) 

_  _  Xg  Z  ®  (  ll) 

m  > 


where  m  =  -+-  v. 

Formula  (11)  cannot  be  used  when  the  plane  of  the  AES  orbit 
coincides  or  almost  coincides  with  the  plane  of  the  earth* s  equator. 
The  method  for  determining  the  AES  orbit  described  in  this  work  Is 
inapplicable  in  these  cases.  As  is  known,  for  Soviet  AES*s  i  =  65°. 

To  determine  the  three  unknowns  _a,  £  and  h,  we  can  use  the  three 

Lambert  equations: 


where 


nx,  =  (i,  —  sin  t.)  —  (l.  —  sin  £<), 


•  ) 


(12) 


=  xt  —  Xj,  r)  =  x)-*-yiJ-*-ziJ> 
byt  =  yt—yj,  rl  =  x\-*-y\-*-z\^ 

—  zt~  zj,  cJ  =  Xx*-*-  Ay]  -+■  Az *. 


k  =  107.0892  (a  value  corresponding  to  the  Qauss  constant  in  the 


% 
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problem  of  the  determination  of  planetary  orbits)  . 


The  subscripts  i,  J_,  and  k  =  1,2  and  3  must  also  satisfy  the 

conditions  k  >  j  and  k  ^  J  ^  i. 

ei 

The  +  or  -  sign  for  sin  — are  assumed  to  be  dependent  on  the 
fact  that  the  AES  makes  an  even  or  odd  number  of  complete  rotations 
about  the  center  of  the  earth  in  the  interval  of  tine  between  the 

6  i 

and  k  observations.  The  choice  of  the  sign  sin  is  made  in 
accordance  with  the  rules  which  were  indicated  in  the  work  of  Subbotin 
[if]  (1941).  The  same  work  also  gives  the  selection  rules  for  that 
quadrant  in  which  it  is  necessary  to  take  the  value  . 

The  system  of  three  equations  (12)  can  be  solved  by  Newton's 
method,  if  the  approximate  values  of  the  unknowns  a,  £  and  h  are 
known. 

The  corrections  to  the  unknowns  Aa,  Ag  and  Ah  satisfy  the  system 
of  linear  equations: 

'  Ai&a Gi$g = 

in  which 

Ff  —  nXf  —  8(  -+-  sin  «< —  sin  8,-, 

A‘~  —  ~  U-5  nrf  —  Mt(rj -*-rk  <»<)-+- /V4(r,-+- r* 

Gi= — —  (15) 

=  —  "jt*  =  Fijyj  ■+-  RuS/t • 

Here 


=  Pi  —  Nt  —  M{, 

Q  =  Ni  +  Mi, 
Fij=-^(P<q/-*-Q<q<j), 

( p>q *  — 


(15a) 
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where 

l  r, i  Li Nt  r  . 


q*j=Tt  (Ax<h Az.v>)- 


In  formulas  (12),  (15).  (1*0  and  (15)  all  the  angles  are  expressed 
In  radians. 

To  determine  the  approximate  values  of  the  unknowns  a,  g  and  h, 
we  recommend  this  method.  If  we  assume  that  the  orbit  of  the  AES 
Is  circular,  then  rx  =  r2  =  r3  =  a  and  Lambert »s  equations  are  ex¬ 
changed  for: 


xjxk'*-viy<c*ifx  _ 

o* 


4t< 


COS  — 7=-  . 
a  Vfl 


(16) 


The  coordinates  x,  j_,  and  £  are  calculated  from  formulas  (6); 
p  appearing  In  these  formulas  Is  determined  from  the  last  of  formulas 
(8)  In  which  r  =  a. 

Equation  ( 16)  may  be  solved  graphically.  It  Is  necessary  to 
construct  a  curve  from  the  parametric  equations: 

-*7= =£.(*)• 
a  v  a 


This  curve  at  high  values  of  kT1  Is  almost  linear  and  the  scale 
of  values  of  parameter  a  on  It  is  almost  uniform.  The  points  of 

Intersection  of  this  curve  with  the  curve  corresponding  to  equation 

t;  =  cos 

are  the  desired  values  of  the  unknown  a.  After  solving  equations  ( 16) 
we  will  have  three  systems  of  values  of  the  unknown  a.  For  the  approx¬ 
imate  value  of  the  major  semiaxis  of  the  AES  orbit  _a  we  take  the  num¬ 
ber  close  to  the  almost  coinciding  solutions  in  the  three  systems  of 
solutions  of  equations  ( 16) . 
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The  approximate  values  of  the  unknowns  g  and  h  are  obtained 

FIRST  **  a  — 

by  the  method  of  least  squares  from  the  system: 

*ig  -*-!Jih  +  zt  =  0  (/ =  1,2,  3). 

The  coordinates  y.^  and  are  determined  from  formulas  (6)] 
pi  is  found  from  the  last  of  formulas  (8)  in  which  for  r^  we  take 
the  above  adopted  approximate  value  of  the  major  semiaxis  of  the  AES 
orbit . 

After  we  have  found  the  approximate  values  of  all  three  unknown 
values  aQ,  gQ  and  hQ  from  equations  (14),  we  determine  the  corrections 
for  these  approximate  values  AaQ,  AgQ  and  AhQ. 

The  new  approximate  values  of  the  unknown  values 
a1  =aa-+- Aa0,  *,  =  &-*-  Ag0,  A,  =  /»„  -+-  AA0 

( first  approximation)  are  substituted  in  equations  ( 12)  and  if  they  do 
not  satisfy  these  equations  with  the  required  precision,  then  using 
equations  ( 14)  we  find  the  corrections  Aat,  Agi  and  Ahi  in  the  second 
approximation  satisfies  equations  (12)  with  the  necessary  precision. 

The  remaining  elements  of  the  AES  orbit  are  found  in  the  follow¬ 
ing  way.  The  set  of  formulas: 


V.  =  1 - 7- ,  ecos Et=qt, 

i _ R.  ninir  —  9jCoi2^i—  <j3 

Ei—E2—2g»,  Ei  —  Et-+-2gl, 

M(  =  E,—e  rin  £.■  =n{ti  —  /.) 


) 


(17) 


enables  us  to  find  the  eccentricity  _e  of  the  AES  orbit  and  the  moment 
the  AES  passes  the  perigee  of  its  orbit  t^. 

The  matrix  of  the  projected  coefficients  are  obtained  from 
formula 


b 

* 
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[Vsl  = 


where 

i  rlWSl  LINE 


OF 


TtX 


The  element  cd  Is  the  angular  distance  of  the  perigee  of  the  AES 
orbit  from  Its  ascending  node  at  the  equator  and  is  determined  from 
equations: 

P,  =  sin/sinw,  Qt  =  sin :  eos  a>.  (19) 

During  the  calculations  it  is  desirable  to  check  the  results  in 
passing. 

Appendix 


Determination  of  the  Orbit  of  the  Rocket  of  the 
Third  AES  (1953  6,)  from  Three  Observations  on 

29  July  1958. 

All  the  calculations  given  below  (1-7)  were  made  by  M.  V.  Kuznetsova. 

1.  Observation  Data  Taken  From  Bulletins  of  the  AES  Visual  Observation 

Station 


1 


2 


3 


Tashkent,  Astri 
nomical  Observ¬ 
atory  of  the 
Ah  Uzbek.  00R 
(Bulletin  ho. 

4,  P.  19) 


■  Pulkovo, 
hiain  Astro¬ 
nomical  Ob¬ 
servatory, 

(Bulletin  ho, 

10,  p.  21) 


Kiev,  hiain 
Astronomi¬ 
cal  Observa¬ 
tory,  Ah  biooK 
(Bulletin  i«Os 

9,  p.  20) 


UT  . 

•iBM.O 

4ISJ0.0 


0.826447 

4686079 

50.6283 


0.S93813 
25080400 
-*-77.081 1 


0.969969 

33484738 

-•-15.1700 


2.  Auxilliary  Values  Depending  on  Basic  Values  Obtained  from  Obser¬ 
vations 


1 

2 

3 

r'  cot  f'  . 

0.752065 

0.504727 

0.638829 

. 

0.656881 

0.860386 

0.766170 

A::::::: 

1.000078 

1.000012 

1.000029 

31387C88 

2998 120S 

32687825 

v . 

0.1098 

0.1098 

0.1098 

S—  V . 

313.6590 

299.0110 

326.6727 

X  . 

0.43579 

—0.07632 

0.87094 

1* . 

0.46096 

—0.21014 

—0.41591 

>  . . 

0.77305 

0.97469 

0.26168 

X . .  . 

0.51924 

0.24478 

0.53379 

Y . 

— 0.54413 

—0.44140 

—035100 

(  cont .) 


* 
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-  —  - 

1 

2 

3 

z  .  .  . 

0.65693 

0.86040 

0.76619 

HR  ST  ft 

0.99862 

0.997S2 

0.997S9 

R*  \ . 

0.99725 

0.99504 

0.99518 

C 

0.48330 

0.91269 

0.81138 

s* . ■ .  s 

0.76368 

0.16203 

0.33684 

in  twenty  four  hour  periods, 

0.0761  S6 

0.143522 

0.067366 

in  radians) . ' 

8.1S549 

15.36966 

7.21417 

Graphic  Determination  of  the  Major  Semiaxis  of  the  Orbit 


i 

!  2 

3 

,0—1.1 

p  . 

0.18477 

0.11101 

0.12305 

X . 

0.59976 

0.23631 

0.64096 

—0.45896 

—0.46473 

—0.40218 

X  . 

0.79977 

0.96890 

0.79839 

0.91875 

0.99796 

0.93362 

e  : . 

7.06905 

13.32218 

6.25313 

p  . 

[  0.33909 

0.21778 

0.23S93 

X  . 

0.66701 

0.22816 

0.74188 

. . 

—0.38782 

—0.46716 

—0.45037 

Z  . . 

0.91906 

1.07267 

0.82871 

?  ::::::::: 

0.88722 

0.99385 

a92150 

6.20411 

11.69213 

S.4S802 

The  drawing  gives  the  points  of  intersection  of  the  cosine  curve 


■q  =  cos  £  and  the  straight  line  connecting  the  points  with  the  coor¬ 
dinates  i  and  T]  calculated  above  for  a  =  1.1  and  a  =  1.2.  It  follows 
from  the  fact  that  the  scale  of  values  of  a  on  the  straight  line  is 
almost  uniform  that  we  may  find  the  values  of  a  at  the  points  of 


intersection.  In  the  given  case  for  all  three  combinations  of  the 


three  observations,  two  at  a  time,  we  obtain  the  approximate  general 
solution  a  =  1.14.  The  general  solution  is  also  taken  to  be  the  zero 


approximation  of  the  major  semiaxis  of  the  orbit. 


4.  The  Determination  of  £  and  h  in  zero  approximation  a  =  1.14 


i 

2 

3 

t  . 

0J4877 

0.1S3S8 

0.16982 

X . 

0.62765 

0.23304 

0.68169 

. . 

—0.42946 

—0.47374 

—0.42163 

*  . 

0^4924 

1.01038 

0^1063 

•  *  —  107.0892 
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FIRST  LINE  OF  TEXT 


0.91295  g0  — 0.66737  h0  =  — 1.32108 
—0.66737  g0  -+-0.58664  A0  =  1.18516 
go= 0.17679 
A0  =  2.22138 


5.  The  Determination  of  the  Corrections  AaQ,  AgQ  AhQ  ^ero 

Approximation  and  Corrections  Aaj,  Agi  and  Ahj  to  the  first  Approxima¬ 
tion 

The  calculation  diagrams  are  the  same  in  both  cases.  Below  we  pre¬ 
sent  the  complete  calculation  for  the  second  case 


1 

2 

3 

m . 

1.84375 

0.50559 

— 0.50801 

p  . 

0.24211 

0.12266 

0.16804 

x . 

0.62475 

0.23542 

068014 

y . 

—0.43253 

—0.46718 

—042089 

X  . 

0.84409 

0.97996 

0.81016 

rJ  .  .  .  :  . 

1.28988 

1.23400 

1.29610 

r  . 

1.13573 

1.11C86 

1.13846 

Ax . 

0.44472 

0.05539 

—038933 

. 

0.04629 

0.01164 

—0.03465 

Ax  ......... 

—016980 

—0  03393 

0.13587 

9  . .  . 

0.47828 

0  06599 

0413S1 

rr4-/-V|  .  ,  .  , 

2.72760 

2.34018 

2.66040 

r  +  r'-. . 

1.77104 

?.2k.  82Q 

1.83278 

s 

»ln  y . 

-0.77333 

a71631 

076374 

•J 

«U-y  .......  . 

0.62314 

i 

a69S81 

—0.63391 

C 

2  •  . . .  • 

4.02567 

7.08168 

2.27251 

1 

2  . 

0.67276 

0.76954 

—0.68660 

6*7  •  •  ; . 

1.21977 

1.02656 

—1.18313 

. : 

0.79676 

0.96878 

—0.81964 

AX . . 

6.69827 

12.62342 

5.92515 

—  •  ........ 

—8.05134 

—14.16336 

—4.54502 

i . . 

134552 

1.53908 

— U7320 

»in  . . . 

0.98057 

0.99966 

—098602 

—•Ini  ....... 

—0.97473 

—0.99949 

0.98054 

F . 

—0.00169 

—0.00069 

-♦-0.00147 

M . 

0.534SS 

0.45016 

— 0.51 88 1 

N . 

034939 

0.42482 

—0.35942 

P . 

-0.18549 

—0.02534 

0.15939 

Q . 

0.88427 

0.87498 

— 0.87823 

0.63871 

0.93204 

0.86030 

—0.43734 

0.04962 

—0.19478 

0.67667 

0.52313 

040943 

::::::: : 

—1.10685 

0.01477 

014799 

**  . 

1.49197 

0.94393 

100503 

A  ~ . . 

8.07493 

16.50528 

8.42747 

C . 

0.75417 

065123 

0.32906 

H . 

—0.11086 

—0.40368 

— 0.533 S4 
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The  corrections  to  the  first  approximation: 

=  —0.000099, 

=  —0.00200, 

,  Ah,  =  —0.00554. 

»  The  second  approximation  is  taken  as  final: 

a,  =  1.140128, 
ffi  =  0.15334,  . 

A,  =  2.17036. 

6.  Determination  of  the  Values  p,  x,  y,  z,  r,  o,  e,  6  and  F 

With  the  obtained  values  of  ag,  g2  and  h^  we  calculate  the  values 
P>  x,  y,  z,  r,  o,  e,6  and  P  from  a  diagram  identical  with  the  one 
given  above.  In  the  case  being  considered  we  obtained  these  values 
of  the  basic  magnitudes: 


1 

2 

3 

X . 

0.62447 

0.23574 

0.68181 

y . 

— 0i43282 

—046630 

—0.42169 

9  . 

.  084361 

0.97589 

0.81067 

r  . 

1.13533 

1.10696 

1.14012 

0  . 

047777 

0.06706 

0.41198 

nx . 

6.69915 

12.62507 

5.92592 

— 1  . 

—8.05021 

—14.16445 

—4.54760 

1  . 

1.34482 

1.S3928 

— 1.37216 

*in « . 

098381 

0.99963 

— 0.98645 

— *ln  1 . 

—097458 

—0.99950 

0.98334 

F . 

% 

—0.00001 

0.00003 

0.00005 

To  check  the  calculations  we  may  use  the  relation: 

*»  —  s»  =  (*i  —  si )-+-(* j  —  Ss);  (2 y, = 2*r,  -+-  2 gj>. 

The  smallness  of  the  values  F  indicates  that  the  desired  magni¬ 
tudes  obtained  in  the  second  approximation  need  not  be  improved  in 
the  future . 


7.  Elements  of  the  Orbit 

t 

In  conclusion  from  formulas  ( 17)  ,  ( 18)  and  ( 19)  we  obtain  the 
*  remaining  elements  of  the  orbit  and  the  matrix  of  the  projected 
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coefficients . 

HR  ST  USt  Of-  1  Xi 

For  the  1958  5i  satellite  for  the  moment  of  the  second  observa¬ 
tion  i.e.  for  1958,  July  29.893813  U.  T.,  we  obtained  the  following 
elements  of  orbit; 

a  =  1.140128,  e  =  0.07098,  Mt=6. 2  0943, 

/= 659316,  2  =  175.958,  o>  =  34:440. 

and  thi3  matrix  of  projected  coefficients; 

(P.  Q,  R,\  /— 0.83930  0.539S9  0.06404N 

P,  Q»  R,  =  —0.17747  —0.38341  0.90636 
\P.  Q.  R.J  V  0.51388  0.74934  0.41761/ 


The  example  given  above  of  the  determination  of  the  orbit  of  an 
AES  indicates  that  the  method  proposed  by  the  author  makes  it  relative¬ 
ly  easy  to  find  the  AES  orbit  from  three  essentially  arbitrary  obser¬ 
vations. 
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